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Why did Supernova 1054 shine at late times? 
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Abstract. The Crab nebula is the remnant of supernova 
1054 (SN 1054). The progenitor of this supernova has, 
based on nucleosynthesis arguments, been modeled as an 
8—10 Mq star. Here we point out that the observations of 
the late light curve of SN 1054, from the historical records, 
are not compatible with the standard scenario, in which 
the late time emission is powered by the radioactive decay 
of small amounts of ^^Ni. Based on model calculations we 
quantify this discrepancy. The rather large mass of ^^Ni 
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needed to power the late time emission, 0.06 
seems inconsistent with abundances in the Crab nebula. 
The late light curve may well have been powered by the 
pulsar, which would make SN 1054 unique in this respect. 
Alternatively, the late light curve could have been powered 
by circumstellar interaction, in accordance with scenarios 
in which 8 — 10 Mq stars are progenitors to 'dense wind' 
supernovae. 
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1. Introduction 

The Crab nebula is among the best studied objects in the 
sky. Still, the nature of the progenitor star and many as- 
pects of the explosion remain unclear. The ancient obser- 
vations of SN 1054, conducted by astronomers in China 
and Japan, have been analyzed by Clark & Stephenson 
(1977). They conclude that SN 1054 was observable in 
daytime for 23 days and during night for some 650 days 
past explosion. Great attention has been directed to un- 
derstand this light curve in terms of supernova theory 
(Minkowski 1971; Clark & Stephenson 1977; Chevalier 
1977; Pskovskii 1977; Wheeler 1978; CoUins, Clapsy, & 
Martin 1999, and references therein). 

It is now widely accepted that SN 1054 was a core- 
collapse supernova, primarily due to the presence of the 
pulsar. That SN 1054 would have been a Type la super- 
nova (SN la), an old suggestion that was recently aired 
by Collins et al. (1999), is clearly in conflict also with the 



large amount of mass in the filaments, and in particu- 
lar with its hydrogen-rich composition. However, the idea 
that SN 1054 was a normal Type H supernova (SN H), 
as was first suggested by Chevalier (1977), is not unprob- 
lematic. The Crab filaments contain only 4.6 ±1.8 Mq 
of material (Fesen, ShuU, & Hurford 1997) and cruise at 
merely ~ 1400 km s^^ (Woltjer 1972), which give a kinetic 
energy an order of magnitude less than the canonical value 
for core-collapse supernovae, 10^^ ergs. Some of this en- 
ergy may also originate from the pulsar (Chevalier 1977). 
As the early observations of SN 1054 indicate a rather 
luminous explosion, Chevalier (1977) suggested that the 
missing mass and energy of the Crab resides in a hith- 
erto undetected outer shell. That some material indeed 
exists outside the visible filaments was recently shown us- 
ing HST observations (SoUerman et al. 2000). 

Nomoto et al. (1982) constructed a detailed model for 
the progenitor of SN 1054. They argued that the progen- 
itor must have been more massive than 8 Mq in order to 
leave a neutron star and less massive than ^ 10 Mq to 
be consistent with the observed metal abundances. An 
8 — 10 Mq star would eject very little heavy elements 
(Nomoto et al. 1982). In particular, it would eject very 
small amounts of radioactive ^^Ni, responsible for the lat e 
time emission of the supernova ( Mayle fc Wilson 1988 ). 



Correspondence to: jsollerm@eso.org 



In this report we want to draw the attention to the long 
duration of the light curve of SN 1054, and point out that 
this cannot be explained within the standard supernova 
scenario, in which the powering of the emission at these 
late phases is due to radioactive decay of a very low mass 
of s^Ni. 



2. Discussion 

2.1. Supernova Light Curves 

The early light curves of SNe II are powered by the explo- 
sion energy slowly diffusing out of the ejecta. The large 
diversity of light curve shapes in this phase (Patat et al. 
1994) is largely due to variations in progenitor radius, 
envelope mass and composition, as well as in the explo- 
sion energy itself. At later phases ( > 150 days), the light 
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curves of SNe II are often powered by the radioactive de- 
cay of ^^Co— >^^Fe, and then the late evolution is quite uni- 
form (Patat et al. 1994). The ^^Co is itself the decay prod- 
uct of ^^Ni, synthesized in the supernova explosion. The 
late light curve of the well studied SN 1987A was repro- 
duced by models with 0.07 M© of ^^Ni (Kozma & Frans- 
son 1998a, b), and most SNe II hght curves do follow the 
decay rate of ^^Co (Barbon et al. 1984; [Patat et al. 19941). 



In fact, the luminosity on the light curve tail can be 
used to determine the mass of ejected ^^Ni, as has been 
done for several supernovae (see SoUerman, Gumming, & 
Lundqvist 1998, and references therein). 

The radioactive energy from the decay of ^^Co is de- 
posited in the supernova ejecta as 7-rays and positrons. 
The luminosity on the light curve tail from the 7-rays is 
given by = 1.26 x lO^^Ky Mm e-'/^^-^ ergs s"\ where 
is the fraction of the 7-rays trapped in the ejecta, Mni 
is the amount of ^^Ni in solar masses, and 111.3 days is the 
e-folding time for the decay of ^^Co. The positrons con- 
assuming 



4.45 X lO^i Mm e-*/i"-3 ergs s" 



tribute 

they are all deposited locally. 

To determine the nickel mass from the observed fil- 
ter light curves, the 7-ray trapping and the bolomet- 
ric evolution have to be taken into account. In Figure 
1 we use SN 1987A to illustrate these concepts. There 
we plot the bolometric light curve of SN 1987A (trian- 
gles) from Bouchet & Danziger (1993). The upper dashed 
line is give n by the expressions abov e for a nickel mass 
of 0.07 Mq (ISuntzeff fc Bouchet 1990| ) assuming fuU trap- 



ping of the 7-rays {F-y=l). This line fits the light curve tail 
at epochs up to ~ 300 days, but later the observed lumi- 
nosity falls below the luminosity expected for full trap- 
ping, due to an increasing leakage of 7-rays. 

The 7-leakage can be illustrated using a simple model 
with a central radioactive source, where the deposition 
F^~{1 — e^"^), and the 7-ray optical depth evolves as 
T = {t/ti)~^ due to the homologous expansion. Here ti 
is the time when r=l. For SN 1987A, ti ~ 610 days gives 
a reasonable agreement with the bolometric magnitudes 
as shown by the dotted curve in Figure 1. This means 
that ^ 60% of the 7-rays were trapped at 650 days for SN 
1987A. 

Also plotted in Figure 1 is the V^-band light curve of SN 



1987A (square symbols, Suntzeff & Bouchet 1990). It can 



be seen that this filter light curve deviates substantially 
from the bolometric light curve after about 500 days. At 
these late phases, the ejecta temperature is low and most 
of the fight is instead emitted in the infrared (Kozma & 
Fransson 1998a, b). 

2.2. The late time luminosity of SN 1054 

The Crab is located at a distance of about 2.0 kpc 
([Trimble 197^) with an extinction of E(B ~ V) ^ 0.52 
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Fig. 1. 

The upper long dashed line labeled '0.07' gives the abso- 
lute bolometric magnitude expected from the radioactive 
decay of 0.07 Mq of ^^Ni, assuming fufi trapping. The 
dotted line allows for 7-ray leakage with ti — 610, as ex- 
plained in the text. This gives a decent fit to the bolomet- 
ric magnitudes of SN 1987A (50 kpc, E{B ~ V) = 0.16) 
shown by the triangles (Bouchet & Danziger 1993). The 
square symbols show the absolute y-band magnitudes for 
SN 1987A (Suntzeff & Bouchet 1990), which deviate sub- 
stantially from the bolometric curve at late phases. The 
lower dashed line is the absolute bolometric magnitude ex- 
pected from the radioactive decay of 0.006 M© ^^Ni, again 
assuming full trapping. The horizontal line represents the 
naked eye detection limit for a source at a distance of 2.0 
kpc and with an extinction of E{B — V) = 0.52, the pa- 
rameters for SN 1054. The error bar shown at day 650 for 
this detection limit encapsulates the errors described in 
the text. 



faded from visibility some 650 days after discovery. The 
detection limit for night time observations was estimated 
to be 5.5 mag by these authors (Clark & Stephenson 1977). 
SN 1054 must thus have had an absolute V (actually vi- 
sual) magnitude of My = —7.6 at this epoch. In Fig- 
ure 1 we have indicated this limit together with an error 
bar that encapsulates distan ces in the range 1500 — 2200 
pc ( Davidson fc Fesen 1985 ), and an error of ±0.04 in 
E{B — V). As a detection limit of 5.5 magnitudes may 



(SoUerman et al. 2000). From the historical records ana- 
lyzed by Clark fc Stephenson (1977) we know that SN 1054 



be regarded too high, we encapsulate it with conservative 
limits of 5.5^q'3. This is also included in the error bar in 
Figure 1. 

The Crab progenitor has been modeled as an 8—10 Mq 
star ( [Nomoto et al. 1982 ). The amount of ^^Ni ejected in 
such an explosion is supposed to be very small. Detailed 
calculations performed by Mayle fc Wilson (1988) indicate 
that no more than 0.002 Mq of ^^Ni should be ejected 
from supernovae with progenitors in this mass range. 
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The ^^Ni eventually decays to ^^Fe and thus the cur- 
rent amount of iron probes the mass of ejected nickel. 
Abundance analyses in the Crab are nontrivial, but in- 
dicate a subsolar iron abundance (Davidson 1978,1979; 
Henry 1984; Bautista et al. 1996). For an ejecta mass of 



4.6 Mq, solar abundance corresponds to 0.006 Mq of iron. 
This is consistent with a low mass of ejected nickel, as sug- 
gested by the explosion models of Mayle & Wilson (1988). 

Here we simply want to point out that the naked eye 
observations of SN 1054 at 650 days seem to be incon- 
sistent with the standard scenario where supernovae from 
8 — 10 Mq stars are powered at late times by the radioac- 
tive decay of very small amounts of ^^Ni. 

Even assuming full trapping and that all of the emit- 
ted flux emerges in the visual band, {Fj—1, Mv=Mhoi), 
just barely keeps the supernova shining at 650 days past 
explosion. The lower dashed line labeled '0.006' in Figure 
1 shows a full trapping, bolometric case for this amount 
of ejected ^^Ni, and is fairly close to the lower limit of the 
luminosity needed for naked eye visibility of the supernova 
at 650 days. 

This discrepancy is of course much more pronounced in 
the more realistic case when the conservative assumptions 
above are relaxed. For example, assuming that SN 1054 
had the same 7-ray leakage and bolometric correction as 
SN 1987A, ~ 0.05 Mq of ^^Ni would be required to reach 
Mv = -7.6. 

2.3. Modeling 

Clearly, it is too simplistic to directly compare SN 1054 
with SN 1987A. To quantify the discrepancy mentioned 
above, we have therefore modeled SN 1054 in some de- 
tail. Lacking an explosion model, we base our input 
model on the available observations. We have adopted 
an ejecta mass of 4.6 Mq (Fescn et al. 1997) and a max- 



imum velocity of 2300 km s"^ (Clark et al. 1983; 



Davidson & Fesen 1985). We assume that the density and 



composition are constant throughout the ejecta. We use 
solar abundances from Cameron (1982) except for helium, 
for which we use a ratio of the mass fractions of helium 
and hydrogen, X(He)/X(H)=2 (e.g., Henry 1986). 

A central source of radioactive ^^Ni is assumed. Cal- 
culations have been done for 0.006 Mq, 0.04 Mq and 
0.07 Mq of s^Ni. The emission from SN 1054 is modeled in 
detail from 200 to 800 days. The code is fully described in 
Kozma & Fransson (1998a, b). In the decay of ^^Co, 7-rays 
and positrons are emitted, and the thermalization of these 
is calculated using the Spencer- Fano formalism (Kozma & 
Fransson 1992). This time-dependent code successfully re- 
produces the detailed late time observations of SN 1987A 
(Kozma 2000). 

With this model, we estimate that only 32% of the 
7-rays were trapped at 650 days for SN 1054. The low 
trapping is of course due to the low mass of the ejecta, 
as compared to SN 1987A. In Figure 2 we plot the light 
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Fig. 2. Absolute V^-band Hght curves for the Crab model. 
The upper curve is for 0.07 Mq of ^^Ni, while the lower 
curve is for 0.006 Mq of ^"^Ni. The detection hmit is the 
same as in Fig. 1. 



curves for the SN 1054 model for two different amounts 
of ejected ^^Ni. At 650 days, a smaller amount of the lu- 
minosity comes out in the V-band for the low nickel-mass 
model. This is because low amounts of nickel means less 
heating, and as the temperature of the ejecta decreases 
the emission is pushed further into the infrared. 

From Figure 2 we directly see that 0.006 Mq of ^^Ni 
could not have provided the observed luminosity at 650 
days. Instead, the model with 0.07 Mq of ^^Ni (similar 
to the mass ejected in SN 1987A), comes close to a V- 
magnitude of 5.5 at day 650. Within the framework of 
this model, a mass of ^^Ni of 0.06t°:°^ Mq was required 
to power the supernova light curve when it faded from 
night time visibility. This includes the errors estimated 
above for distance, reddening and naked eye-sensitivity. 
We have done the same calculations also for a filter sim- 
ilar to the eye sensitivity (Rhodopsin absorption curve, 
Kitchin 1991), but the deviation from the F- filter is small, 
0.07 mag for the model with 0.07 Mq of ^^Ni at 650 days. 
Finally, although according to Clark & Stephenson (1977) 
there is no reason to question the dates of the Chinese 
sightings, an error of ±20 days in duration would intro- 
duce a 30% error in flux. We note that suggestions for 
alternative durations, although often based on less direct 
sources (e.g., Pskovskii 1977; Colhns et al. 1999), usually 
argue for longer visibility of SN 1054, thus enhancing the 
problem of keeping it shining at late times. 

3. Powering scenarios 

3.1. Radioactive decay 

If radioactivity alone powered the late time emission, the 
required mass of ^^Ni estimated above, 0.06tQ Q3 Mq, is 
significantly higher than that obtained in the explosion 
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models by Mayle & Wilson (1988), 0.002 Mq. Although 
the unknown mechanism of core-collapse explosions makes 
the existing explosion models rather uncertain, including 
the exact amount of nickel mass, there are other reasons 
to believe that low mass SNe eject very small amounts of 
^^Ni. Galactic chemical evolution models imply that the 
amount of iron produced by SNe II is already quite high, 
and Samland (1998) even suggested that the lower mass 
limit of stars producing iron should be 11 — 12 M©, instead 
of the conventional 10 M©. 

Nevertheless, the above arguments cannot conclusively 
argue that SN 1054 did not eject the 0.06 M© of 
required to maintain the late time luminosity. We there- 
fore turn the attention to the measured abundances in the 
Crab nebula. 

Actually, only helium is enhanced in this supernova 
remnant, while the abundances of C, N, O and Fe seem 
to be (sub)solar ( pavidson fc Fesen 1985). In §2.2 we saw 



that solar abundance of iron for 4.6 M© ejecta corresponds 
to 0.006 M©. The amount of ^^Ni required to power the 



late time light curve, 0.06_q;o3 M©, would thus correspond 
to a present iron abundance which is 9lg times higher than 
solar abundance. 

If radioactivity powered the supernova, the iron must 
be locked up in dust. Dust is clearly present in the nebula, 
and it is nontrivial to determine the dust mass (Sankrit 
et al. 1998), although estimates seem to indicate that the 
dust mass is low (Davidson & Fesen 1985). Note that this 



scenario would also require large amounts of other met- 
als to be locked up in dust. As the iron-group elements 
are produced in the very center of the exploding star, 
we do not expect to find a large amount of iron, while 
the CNO-abundance has not been enhanced. For example, 
SN 1987A, which ejected 0.07 M© of ^^Ni, also ejected al- 
most 2 solar masses of oxygen. Even the 11 M© model of 
Woosley & Weaver (1995), which ejects 0.07 M© of ^^^Ni, 
and only 0.136 M© of oxygen, provides an oxygen abun- 
dance higher than solar. We believe that the scenario in 
which SN 1054 was powered by 0.06 M© of ^'^Ni, but where 
all the nucleosynthesized metals are presently locked up in 
dust, is not very convincing. We therefore briefly turn our 
attention to alternative scenarios. 

3.2. The Pulsar 

If radioactivity did not power SN 1054 at late times, the 
obvious candidate is the pulsar that powers the nebula 
today. Already Chevalier (1977) suggested that the pulsar 
could contribute to the late time supernova light curve. 
In fact, unless the bulk of the iron is locked up in dust, 
and circumstellar interaction was unimportant (see §3.3), 
the pulsar would have to make up for all of the late time 
luminosity, as 0.006 M© of ^^Ni could not have kept the 
SN visible for more than ^ 500 days (Fig. 2). 

There are many ways in which a pulsar can contribute 
to the optical luminosity of a supernova. Accretion onto 



the neutron star could be either spherical or in a disk, de- 
pending on the angular momentum of the infalling matter. 
A steady spherical accretion is supposed to reach a maxi- 
mum Eddington luminosity of LEdd = 3.5 x lO'^* ergs s~^, 
and Chevalier (1989) suggested that this scenario might 
be responsible for the late time luminosity of SN 1054. He 
worried, however, about the fact that the accretion lumi- 
nosity would not be able to escape from the vicinity of the 
neutron star during the first months, as the luminosity is 
trapped by the inflow (see also Benetti et al. 2000). As 
shown in Figure 2, this need not be a problem, as even 
a low mass of ejected ^^Ni is able to keep the supernova 
shining for some 500 days. If the accretion is mediated via 
a disk, the luminosity could perhaps be even larger than 
in the spherical case. 

Unknown in these scenarios is the fraction of the lu- 
minosity escaping in different bands. In fact, the full Ed- 
dington luminosity would have to fall in the visible band to 
make SN 1054 observable at 650 days (3.5 x 10^* ergs 
corresponds to Mboi = —7.7). This appears to be rather 
unlikely. 

Pulsar nebulae in supernovae, with special attention 
to the Crab, were investigated by Chevalier & Fransson 
(1992). Unfortunately, they did not address the question 
highlighted here, the luminosity of SN 1054 at 650 days. 
In their scenario, the pulsar powers a bubble that shocks 
and ionizes the supernova ejecta, and their calculations 
suggested that 1.5% of the total pulsar luminosity could be 
converted to radiation. This would be too low to account 
for the luminosity of SN 1054 at day 650. However, if 
the pulsar bubble is a significant source of synchrotron 
emission (Chevalier 1996), the efficiency might be higher. 
Another interesting possibility is that the pulsar was born 
with a very rapid spin period (see Atoyan 1999), and thus 
a very high initial spin-down luminosity. 

In summary, although not investigated in detail, there 
are several ways in which a pulsar could contribute to 
the luminosity of supernovae. This has so far never been 
unambiguously observed. If the pulsar was indeed respon- 
sible for the late light curve, this would make SN 1054 a 
unique case. 

3.3. Circumstellar interaction 

Another mechanism important for supernova light curves 
is interaction with circumstellar material (CSM) . Circum- 
stellar interaction is responsible for the emission at late 
phases for a number of supernovae, such as SNe 1979C, 
1980K, 1988Z, 1993J and 1995N. Such interaction can 
maintain a SN luminosity of severalx 10'^^ ergs s^^ for 
many years (e.g., Chugai, Danziger, & Delia Valle 1995). 
In this respect, we note that the Crab could fit into the 
scenario of 8 — 10 M© stars being progenitors to the so 



called 'dense wind' supernovae (Chugai 1997), where the 



supernova ejecta interact with a dense superwind from the 
progenitor star. 
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The dense wind must have extended out to ^ 6 x 10^^ 
cm for a maximum ejecta velocity of Woj ~ 10"* km s~^. 
This means that the wind started > 2{V(,]lv^) years be- 
fore the SN breakout. Here v^j is the wind velocity. Com- 
paring with the line fluxes in Chevalier & Fransson (1994), 
we find that the late light curve can be explained ii M/v^ 
was typically in excess of 10~^ Mq yr~^/10 km s~^. Here 
M is the mass loss rate during the superwind. The swept 
up wind should in this case now coast freely outside the 
present nebula. However, no sign of ejecta, or swept up 
wind, moving at velocities of order 10* km s^^ has been 
identified (e.g., Fesen et al. 1997). 

Chugai & Utrobin (1999) suggested that SN 1054 was 
similar to the low-energy (~ 4 x 10^° erg) event SN 1997D 
(Turatto et al. 1998; Benetti et al. 2000), and that no fast 
ejecta exist. In this case, circumstellar powering of the late 
light curve could be more difficult as the luminosity of the 
circumstellar shock scales as oc v^-, . Furthermore, in a low- 
energy explosion, the peak luminosity is more likely to be 
lower, and Chugai & Utrobin (1999) therefore suggested 
that circumstellar interaction was important for SN 1054 
also at the early phase. A low-energy explosion model with 
a circumstellar shell was explored already by Falk & Ar- 
nett (1977), and they found a high efficiency in the con- 
version of shock energy to visual light during the peak, 
in accordance with the suggestion of Chugai & Utrobin 
(1999). Although there is thus no direct support for cir- 
cumstq llar interaction as the cause for the late emission in 
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the Cr ab, it is still premature to rule out this possibility. 



4. Summary 

SN 1054, the creation of the Crab nebula and the Crab 
pulsar, is the typical example of an 8 — 10 Mq super- 
nova. Such supernovae are expected to eject only minute 
amounts of ^^Ni, and here we demonstrate that SN 1054 
could not have been powered by such a small mass of nickel 
at late times. The required amount of ^^Ni, ~ 0.06 Mq, 
is much larger than suggested by abundance analyses of 
the Crab, and we therefore discuss alternative solutions. 
The pulsar may have powered the supernova, which would 
make SN 1054 unique in this respect. Alternatively, the 
progenitor of SN 1054 could have had a dense wind and 
the supernova could then have been powered by circum- 
stellar interaction. 
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